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We examine nonlocal conductivity in high-temperature superconductors from a phenomenological 
point of view. One wants to deduce the properties of the conductivity, especially its inherent length 
scales, from the transport data. Although this is a challenging inverse problem, complicated further 
by the experimental data not being completely self-consistent, we have made some progress. We 
find that if a certain form for the conductivity is postulated then one requires positive "viscosity" 
coefficients to reproduce some of the experimental results. We are able to show that the effects of 
surfaces on the conductivity are likely to be important and draw comparisons with the treatment of 
the surface within the hydrodynamic approach put forth by Huse and Majumdar. We also develop 
an approximation scheme for the conductivity which is more robust than the hydrodynamic one, 
since it is stable for both positive and negative viscosity coefficients, and discuss the results obtained 
using it. 

PACS numbers: 74.60.-w, 74.25.Fy, 74.60. Ge, 74.20.De 



I. INTRODUCTION 



The measurement of a substantial nonlocal conductiv- 
ity in high-temperature superconductors in a magnetic 
field is thought to imply the existence of moving vortex 
lines having coherence lengths of the order of the sample 
thickness, as opposed to pancake vortices readily slid- 
ing past one another (for a general review, see Blatter 
et al.u). Measurements of nonlocal effects probe the in- 
herent length scales of the problem and thus can be used 
to investigate issues such as whether the-.decoupling and 
melting transitions occur simultaneouslju. The claims of 
Safar et al. H to have observed a sizable nonlocal effect in 
twinned YBa2Cu307_,5 (YBCO) are based on two sets 
of measurements. In the first, which we refer to as the 
top geometry, a current is put into and drawn out of the 
top of a modified flux transformer while the voltage dif- 
ferences, V top and Vbot, are measured (see Fig. |l|). In the 
second, the side geometry, the current is withdrawn from 
the bottom and tha-voltages, Vieft and V r i g ht, are mea- 
sured. Safar et al.u find that the ratios Vbot /Vtop and 
V r ight /Vie ft both approach one as they near the melting 
transition. Taken individually either result might be ex- 
plained by a local though anisotropic conductivity; but 
taken together the results are inconsistent with a local 
description. Safar et al.u confirm this by analyzing each 
data set as though the conductivity were local (the Mont- 
gomery analysisu) and extracting from each the apparent 

conductivity ratio ai a J / 1 erf 1 ) ', finding a huge discrepancy 
in these apparent ratios. 



Iin f A I out iin 

t y^x k t ^ T 4 



B 





fG. 1. The modified flux-transformer set-up of Safar et 
The diagram on the left shows the arrangement of the 
terminals for the top geometry, and the right shows the side 
geometry. The samples are single crystals of YBCO with 
the magnetic field aligned along the c-axis of the crystal. 
Vi e ft = V% - Ve and V right = V 3 - V 7 . 



Eltsev and Rapp@ dispute the Safar et al. claim. 
They performed similar measurements but did not see 
V r ight /Vieft — * 1- On the other hand, they may have 
seen nonlocal effects in a tilted geometry in which the 
current is extracted from terminal 6 (instead of terminal 
5) and the ratio (V3 — Vs)/(V~2 — V7) measured. In a co 
parison of twinned and untwinncd YBCO, Lopez et 
find that in the untwinned YBCO, the strongest signa- 
ture of nonlocality seen by Safar et al. □, Vbot ~ Vtop-, is no 
longer found for any significant temperature range above 
the "melting" transition. A feature of our studies below 
is that substantial nonlocal effects are only present when 
a characteristic length (presumably the phase coherence 
length) is of the order of the sample thickness. Now the 
phase coherence length scale along the field direction (ac- 
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cording to Ref.Q) grows exponentially rapidly as the tem- 
perature is lowered in such a way that the temperature 
interval over which nonlocal effects might be visible in 
the vortex-liquid region is only perhaps within 0.3K of 
the temperature at which pinning drives the resistance 
rapidly to zero. This is of the same order of magnitude 
as the rounding of the zero-field transition due to sample 
inhomogeneities and as a consequence it will be hard to 
disentangle the various effects Jrom each other in the un- 
twinned results of Lopez et alB. If one supposes that the 
long length scale causing the nonlocality in thcptwinncd 
case is caused by a Bose-glass-like mechanisnJffl, then 
one would expect the coherence length to increase only 
as a power law, and the width of the temperature interval 
over which nonlocal effects arc visible may therefore be 
wider in the twinned case. 

The assertions of nonlocal effects in Bi2Sr2CaCu20 x 
(BSCCO) are less dramatic than those in YBCO. In theix 
measurements on single-crystal BSCCO, Keener et al. E£l 
never observe the ratios Vbot/Vt op and V r i g ht/Viejt simul- 
taneously approaching one. Nevertheless, when they per- 
form a Montgomery analysis on their data, they do see 
discrepancies in the apparent p.tio a x a x /crfz ■ Conversely, 
measurements by Buach et al. Ell on single-crystal BSCCO 
and by Doyle et al. 113 on BSCCO with columnar defects 
are claimed to be consistent with local resistivity. These 
seemingly contradictory results could be, caused by ap- 
proximations used in the local analysisliil and might be 
resolved by using the full analysis or better suited, ap- 
proximations, such as the one proposed by LevirJlj. A 
theoretical framework which could provide some quanti- 
tative analysis of these results — for example, by the ex- 
traction of a temperature-dependent length scale — would 
obviously be helpful in the interpretation of these and 
other results. Our aim is to develop such a framework. 
We will approach the problem phenomenologically, at- 
tempting to relate the current-voltage characteristics to 
the form of the (nonlocal) conductivity. 

When a material has a nonlocal conductivity, the ap- 
propriate form of Ohm's law is given by 

J»= / ^(r,r')£,(r')dV, (1) 

where symbols have their usual meaning in this context, 
and the integral is taken over the volume of the sample. 
For a nonlocal conductivity, oy^r, r') ^ 8(v — r'). In 
momentum space, for a translationally invariant system, 
this relation becomes 

j M (k) = ^(k)4(k), (2) 

where j M (k) is the Fourier transform of j M (r), and sim- 
ilarly for the other quantities. It should be noted that 
nonlocal effects will only be observable when the length 
scale of the nonlocality is of the same order or larger than 
the distance between leads. 

The best known theoretical work on the subject is the 
"hydrodynamic" approach, expounded upon in general 



by Marchetti et al. E-3 and applied specifically to the con- 
ductivity by Huse and MajumdarEj. This theory is so- 
called because the nonlocal conductivity cr^k) is ex- 
panded in a Taylor series in k, and the expansion is ter- 
minated at order k 2 ; in other words, the conductivity is 
taken to be of the form 

<V/(k) = <v,(0) + r7 M/ 3 7 „ kpkj. (3) 

We revert to the notation of Huse and Majumdar to 
facilitate rGpmparison with that workEd; note that later 
works replace the 77's with S's to prevent confusion 
with the viscosity tensor of the vortex-line liquid, which 
is related but distinctta. Unfortunately, this form for the 
conductivity is unphysical if certain coefficients become 
negative, as shown by Blum and MooreO. Huse and Ma- 
jumdar always assume that the coefficients they use are 
positive. When and whether the coefficients are in fact 
positive or negative will be discussed in more detail in 
Section [|. 

The hydrodynamic analysis leads to a fourth-order par- 
tial differential equation which reduces to Laplace's equa- 
tion in the local limit (r) = 0). It also supplies sufficient 
boundary conditions to solve for the potential V(r). Huse 
and Majumdar argue that there are discontinuities in the 
first derivative of E(r) at the surface. When the conduc- 
tivity (which in hydrodynamics is a differential operator) 
is applied, the result is (5-functions in the current distri- 
bution at the surface, i.e. surface currents. One then 
uses charge conservation, V • j = 0, to translate this out- 
come into boundary conditions on V(r). To handle the 
(5-function it is convenient to integrate over the—surface 
as in the standard Gaussian pillbox argumentaiij — only 
here, because of the surface current, the side surfaces of 
the pillbox contribute even as the volume of the box is 
shrunk down to zero. This gives what initially looks like 
an extra term in their boundary conditions. 

Huse and Majumdar study a two-dimensional geom- 
etry modeling the flux transformer used in the experi- 
ments, the z-axis of which coincides with the c-axis of the 
superconductor. They have performed a detailed analy- 
sis of the situation with one non-zero viscosity coefficient, 
Vxzzx ; which embodies the interaction of pancake vortices 
moving in different ab planes and at different velocities. 
Somewhat surprisingly, their equation and boundary con- 
ditions are symmetric under r) xzzx «-> rj zxxz despite the 
fact that these coefficients would appear to represent very 
different physics. 

As an alterative to the hydrodynamic truncation of 
(^(k), we consider an analysis based on Pade approxi- 
mations to a^v (k) . It incorporates a more realistic large- 
k behavior than the hydrodynamic approach and remains 
solvable. In principle, one can approximate cr^^k) to any 
desired degree of accuracy by using a sufficiently large- 
order Pade approximation. As the order of the approxi- 
mation is increased, our technique of solution continues 
to work, but the computing effort increases rapidly. One 
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stage of the solution involves a partial differential equa- 
tion rather reminiscent of that occurring in hydrodynam- 
ics. In fact, in one instance we can recover the results of 
Huse and Majumdar by means of a limiting procedure 
on the Pade result. 

In the remainder of this paper we first discuss our moti- 
vation for wanting to improve upon and extend the work 
of Huse and Majumdar; this involves an examination of 
whether the relevant coefficients rj in the small-fc expan- 
sion of the conductivity are positive or negative. We 
investigate the current- voltage characteristics in a partic- 
ular geometry (the infinite-slab geometry), which allows 
us to comment on whether positive or negative viscos- 
ity coefficients are needed to explain experimental data 
like that of Safar et al.U (Sec. ||l|). The section following 
that contains details of work using Pade approximations 
to the conductivity. In Sec. [v] we discuss the role of sur- 
faces in determining the conductivity and how the anal- 
ysis of Huse and Majumdar takes account of surfaces. 
Appendices |X| and contain some calculational details, 
while in Appendix y we outline the Bose-glass scaling of 
the conductivities used in some of the numerical work. 



II. ARE THE VISCOSITY COEFFICIENTS 
NEGATIVE OR POSITIVE? 

For stability, the conductivity tensor <7^„(k) must be 
a positive definite matrix. For this to be true of the 
hydrodynamic form, Eq. (|^), certain of the viscosity co- 
efficients rj^apv must be positive; in particular, r\ zzzz and, 
Vxxxx- The work of Mou et a/.ll3 and Blum and MooretZl 
shows that for high temperatures, these coefficients are 
actually negative. Both works use the time-dependent 
Ginzburg-Landau equation as a starting point, so "high 
temperatures" in this context means near the H C 2(T) 
line. Thus to treat nonlocal conductivities in this re- 
gion of the H-T plane, one requires a model that can 
handle these so-called negative viscosities, ruling out the 
hydrodynamics approach. However, we expect substan- 
tial nonlocal behavior occurs only near the melting line, 
where some of the viscosities may very well be positive 
and hydrodynamics a viable approach. 

So what happens to cr^^k) as the temperature is low- 
ered? The arguments of Mou &-,al. Il3 suggest and the 
simulations of Wortis and HuseE£l bear out that as the 
temperature is decreased, rj xxxx changes sign, becom- 
ing positive. Imagine a plot of <J xx {k x ) (see Wortis and 
HusaL3): at high temperatures, a xx (k x ) increases mono- 
tonically els kx is decreased, but at low temperatures, it 
develops a local maximum at a nonzero k x and then falls 
to a finite value (the flux-flow valuccj) at k x = 0. 

In the presence of a magnetic field, the afo-plane con- 
ductivity of a type-// superconductor does not diverge, 
because a current causes the vortices to move, leading 
to dissipation, that is, a non-zero resistance. The con- 
ductivity is thus enhanced if this movement of vortices 



is impeded, for instance by pinning centers. The interac- 
tion of a vortex with other vortices may also inhibit its 
motion. However, for a uniform current the vortices all 
move together and thus their mutual interactions play 
little role in hindering the center-of-mass motion. For 
a nonuniform current, on the contrary, the vortices are 
impelled to change their relative positions, and so their 
interactions do inhibit this sort of motion. Therefore, 
the conductivity may be higher for a nonuniform cur- 
rent. These arguments suggest that at low temperatures 
where the interactions become important the conductiv- 
ity might be higher at some nonzero k x . 

The situation is different for conductivity along the 
c-axis. Like & xx (k x ), at high temperatures, <r zz (fc z ) in- 
creases monotonically as k z is lowered. However, the cur- 
rent is now along the axis of the vortices, the vortices are 
not forced to move, and this time the uniform conduc- 
tivity a zz (k = 0) diverges as the temperature is lowered. 
Hence, there is no compelling reason to expect that a 
local maximum will develop as temperature is lowered, 
and one might suspect that r] zzzz remains negative and 
therefore unsuited for the hydrodynamic prescription at 
all temperatures. This expectation for rj zzzz seems to be 
borne out by low-temperature calculations which have 
the Abrikosov lattice as a starting. pointEil and also by 
the preliminary simulation resultscj. 

The proposition that some of the ry's may be negative 
at all temperatures is one motive for wanting an alterna- 
tive to hydrodynamics; there are others. Hydrodynamics 
is a simple approximation to the actual conductivities, 
but we do not know how good an approximation it is. 
Plus, there is no obvious way to improve upon it — one 
might include terms of order fc 4 , but this necessitates ad- 
ditional boundary conditions to specify the solution, and 
it is unclear what they would be. In the Fourier represen- 
tation of a function, the small-fc terms model well its bulk 
properties, but higher-fc terms are needed to capture the 
behavior at the boundaries. So one might think a proce- 
dure focusing on small-fc would do well in the bulk and 
perhaps less satisfactorily at the surface. But hydrody- 
namics involves a differential equation, and its solution, 
even in the bulk, is determined by boundary conditions, 
i.e. the surface. Thus it is crucial to treat the surface 
properly — even more so, since in the experiments at is- 
sue here, all of the measurements are taken at the surface. 
Huse and Majumdar do take some account of surfaces, 
and in doing so find surface currents, but their theory 
is unable to make any prediction about the length scale 
over which these currents might flow. Such restrictions as 
these provide the impetus to go beyond hydrodynamics. 

There is one further comment that is useful to make 
before proceeding to some concrete calculations. While 
we require that a^ v (k) be a positive definite matrix, this 
does not imply that a^ v {r) must always be positive. Al- 
though conductivities taking on negative values seems a 
littlej-odd at first sight, the simulations of Wortis and 
HusaL3 find that a xx (x,k y = 0) can be negative over a 
range of a few inter- vortex spacings. In fact, any non- 



3 



locality in the conductivity implies that either the real- 
space conductivity or resistivity (or both) take on nega- 
tive values at some points. By definition, the conductiv- 
ity matrix is the inverse of the resistivity matrix, which 
implies that in Fourier space, 6r /xv (k)p J/Q (k) = 8^ a (<5 MQ 
is the Kronecker delta function and the summation con- 
vention is used). For example, if a xy = a xz = we have 
in real space 
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dV a xx (r - r') p xx (r' - r") = 5(r - v"). 



(4) 



-L/2 



L/2 



If r 7^ r", the right hand side of Eq. (Q) is zero; for the 
left hand side to be zero, some cancellation is needed. 
However, there would be no cancellation if cr xx (r — r') 
and p xx (r' — r") are both everywhere positive. Thus, one 
or both has negative regions. For experimental setups 
in which the current is distributed throughout the sam- 
ple, this feature may not manifest rifcsclf in the voltage 
distribution; thus, Wortis and Huseta have proposed ge- 
ometries with very localized currents in order to look for 
it. 



III. THE INFINITE-SLAB GEOMETRY AND 
POSITIVE VISCOSITIES 



The aim of this section is to show that for the choice 
of conductivity given below, it is necessary to have pos- 
itive viscosity coefficients to produce the strongly non- 
local behavior seen by Safar et al.u: Vb ot /V t0 p — ► 1 and 
Vright/Vieft — * 1 simultaneously. Let us consider the two- 
dimensional geometry shown in Fig. |^, with the lateral 
dimension L — > oo; taking this limit eliminates one set of 
boundary effects. Furthermore, let us postulate conduc- 
tivities of the form 



0xx(k) — 5j k ; 
o-* z (k) = a zz (k x ), 



(5a) 
(5b) 



that is, a xx is a constant (i.e. local) and a zz is a function 
of k x alone. The function a zz (k x ) may include a constant 
piece; moreover, that constant, as well as Yt x , may repre- 
sent both superconducting and normal contributions to 
the local conductivity. Restricting the wavevector de- 
pendence to k x alone enables us to solve for the potential 
V(x, z) via Fourier transformation. 



FIG. 2. The geometry used for our calculations of the po- 
tential V(x,z). The c-axis of the superconductor is aligned 
along the 2-axis (called the y-axis in the notation of Huse and 
Majumdar). 

To determine V(x,z) for a given input current, first 
relate V(x, z) to the components of the current 

j x (x, z) = -Sj,. d x V{x, z); (6a) 

p OO 

j z (x,z) = - a zz (x - x')d z V(x',z)dx', (6b) 



using — —d^V. Next note that in the steady state, 
the continuity equation is V • j = 0, which in this case is 

/>OC 

E x dlV(x,z)+ / a zz (x-x')d 2 z V(x',z)dx' = 0. (7) 



Fourier transforming Eq. (^) with respect to x yields 

- k 2 x V(k x ,z) + a zz (k x ) d 2 z V(k x ,z) = 0, (8) 

where we have used the fact that the transform of the 
convolution is the product of the transforms. The solu- 
tion of this differential equation is 

V(k x ,z) = A{k x ) cosh[n{k x )z] + B(k x ) wah[K(k x )z], (9) 
where 



2 s z k 2 

o- zz (k x ) 



(10) 



One determines the functions A(k x ) and B(k x ) from the 
boundary conditions; toward this end, it is convenient to 
Fourier transform the expression for j z (Eq. (pq)) 



d- zz (k x )d z V(k Xl z) = -j z (k x ,z). 



(11) 



Imposing the boundary conditions appropriate for the 
top geometry, namely j z (x, 0) = and j z (x, D) = Jt(x), 
finding V(k x , z) and taking the inverse transform yields 



Vt(x, z) = - 



.Jk, T X 



dk x Jr{k x ) cosh(ftz) 
27r na zz (k x ) sinh(KD) 



(12) 
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Note that charge conservation implies Jr(0) = 0, and 
hence the integral above converges at k x = 0. For the 
side geometry, similar manipulations using the boundary 
conditions j z (x,0) = j z {x,D) = J s (x) give 



V s (x, z) 



dk x Js{k x ) sinh [n ( 



2tt Ka zz (k x ) cosh(^) 



(13) 
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Analytically the rc-axis decay length is controlled by 
the pole structure of the above integrals (see Appendix |A| 
for more details). Rather generically, this length grows if 
the viscosity coefficient (the coefficient of k 2 in the small- 
k x expansion of (J zz {k x )) is positive, leading to features 
such as V r ight/Vi e ft — > 1. Increasing this coefficient (rj) 
also changes the ratio Vbot/Vtop though not necessarily 
in a monotonic fashion. However, there is another way 
to ensure that Vbot/Vt op — ► 1: this is simply to make a zz 
very large, which applies even in the local limit. Thus, to 
obtain results like those of Safar et al.B, we expect that 
a & zz (k x ) which grows large and has a positive viscosity 
coefficient at low temperatures is required. 

These arguments have been checked numerically for a 
variety of conductivities. One example is shown in Fig. |^. 
We have put all of the temperature dependence of the 
conductivity into the length scale £ which is presumed to 
increase as temperature decreases and choose a conduc- 
tivity of the form a zz {k x ) = a z n) + C£ 2 {1 + 2k 2 J 2 )/{l + 
k 2 £ 2 ), which meets the above criterion as I increases. 

(n) 

Note that we have also included a z , a local term that 
does not scale with I, Fig. || shows a plot of Vb t/V t0 p 
and V r ight/Vi e ft as a function of l/£, for this choice, ft 
can be seen that as I — > oo, the two ratios do indeed 
approach one. 

Using conductivities with negative viscosity gives re- 
sults such as those shown in Figs. [I] and |^. In this 
particular case, the conductivity used was a zz (k x ) = 
a z n ^ + C£e~ k * e2 ; however, the results are typical of con- 
ductivities with negative viscosity. Finding negative volt- 
ages is not necessarily unphysical (see the discussion sur- 
rounding Eq. (|J)) but the results clearly do not give us 
the strongly nonlocal behavior seen by Safar et al. a. 



FIG. 3. The ratios of Vbot/Vtop and V r i g ht/Vi e ft for the 
conductivity with positive viscosity coefficients, the form 
of which is provided in the text. The boundary condi- 
tions used were Jt(x) = Jo [S(x — 2L C ) — S(x + 2L C )] and 
Js(x) = —Jq8(x + 2L C ) (4Z/ C is defined to be the distance 
between the current inputs) . The voltages Vt op and Vbot were 
measured at x = ±L C as were the voltages Va, V3, Ve and V7. 
The parameters used are = 5, ai"' = 1, C = 1, D = 1 
and L c = 1; this corresponds to measuring lengths in terms 
of the thickness of the sample and conductivities in terms of 
cri™' , as will also be done in all the other figures. 




FIG. 4. The voltages Vtop and Vbot and their ratio for the 
conductivity a zz (k x ) = ai"^ + Cle~ k:cl , with the same val- 
ues of parameters as in Fig. [| By definition, Vtl is defined 
to be Vtop in the local case for this geometry, which can be 
calculated from Eq. (|A2|) . Note that Vbot/Vtop does not tend 
to 1 as I — > 00. 
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FIG. 5. The voltages Vieft and Vright and their ratio for 
the conductivity a zz (k x ) = u z + C£er k * 1 , with the same 
parameters as for Figs, q a nd H . Here, Vll is defined to be 
Vieft for the local case, Eq. (A2). The notable feature here is 
that V r ight becomes negative. 

Now let us consider the difference between the experi- 
mentaLdata taken from twinned YBCOH and untwinned 
YBCOQ. Recall that for the twinned YBCO, the V top and 
Vbot curves meet at some temperature T t h, and they con- 
tinue toward zero together as the temperature is lowered. 
For the untwinned YBCO, on the other hand, the curves 
only meet just before dropping sharply to zero. We can 
reproduce some of these features by choosing the scaling 
forms for the conductivity appropriately, depending upon 
whether the sample is twinned or untwinned. For the 
twinned case, we use the Bose-glass scaling forms, which 
are discussed in Appendix Although these are sup- 
posed to be valid in the presence of columnar defects, and 
not strictly twin planes, we use them here since there is 
currently no better altetpative, and there is at least some 
experimental evidenced to support a Bose-glass transi- 
tion in twinned YBCO. For the untwinned case, we will 
use the same forms used to generate Fig. ||. 

We choose conductivities given by 



» 



Untwinned: 



Twinned: 



i + kip 



ai n) + C 2 £ 4 



(14) 



& xx = Ox ■ -|- 02* 

(n) di 6 (i + 2kie) 
azz ~ az + TTW 2 ' 



Both choices have the same overall form with positive 
viscosity coefficients, the only difference is in how the 
constants scale with £. 



FIG. 6. The results for Vt op and Vbot as a function of £, us- 
ing the conductivities of Eq. ( Jl4[ ) with parameters ai™' = 5, 
cri™' = 1, D = 1, Ci = 1 and C2 = 5. The upper curve in 
each plot is Vtop, and the lower is Vbot- 

The results obtained are shown in Fig. Firstly, it 
can be seen that in both instances, Vt p/Vb t — > 1 as 
the length scale £ is increased. However, for the twinned 
case, both V top and V uot tend to zero together, whereas in 
the untwinned case, they tend to some finite value. If one 
identifies the sharp drop in the corresponding experimen- 
tal results as the transition from a vortex liquid to solid, 
then this is the correct behavior. The voltages in Fig. ^ 
level off as I grows small, whereas the voltages measured 
in the experiments continue to rise as the temperature 
is increased. This discrepancy is due to our neglect of, 
among other things, the temperature dependence of the 
normal component, which can be nontrivialEJ. Never- 
theless, the results show that varying the scaling behav- 
ior can account for some of the differences between the 
twinned and untwinned YBCO. 



IV. WORK WITH PADE APPROXIMATIONS 

In this section, we lift some of the restrictions in the 
previous section by considering a flux-transformer geom- 
etry with finite lateral dimension L and by generalizing 
the nonlocal conductivities of Eq. (||) to allow either a xx 
or a zz to be nonlocal and depend on either k x or k z . 
As opposed to the infinite-slab geometry, we could not 
make progress with general conductivities, thus we elect 
to use Pade approximations to the conductivities. Two 
important advantages to using Pade forms are: 1) they 
capture the high-k behavior as well as the low-k behav- 
ior of the true conductivity and 2) the resulting equations 
are analytically soluble. In addition, the corresponding 
real-space conductivities can be chosen to be decaying ex- 
ponentials, which, according to Wortis and Huseta, could 
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be the correct form for the high-temperature regime. 
We will consider four different cases: 



rr (\c\ 


— v 


<7xx(k) 


= Z X 




= X X 


o- Z z(k) 





1 + k 2 J 2 ' 

A x 
1 + k 2 £ 2 ' 

1 + k 2 £ 2 ' 



(15b) 
(15c) 
(15d) 



Note that the momentum dependence is upon k z in Cases 
1 and 2 whereas it is upon k x in Cases 3 and 4 and that 
a xx is local in Cases 1 and 4 while a zz is local is Cases 2 
and 3. In the Pade form, £ is the conductivity as k — > oo, 
£ + A is that at k = 0, and I is a length scale. We set 
out to solve V • j = subject to the usual boundary 
conditions on the current at the surface. Each of these 
problems involves an integro-differential equation, which 
can be converted into a partial differential equation with 
linear coefficients, which can in turn be solved by sep- 
aration of variables. This procedure is outlined in Ap- 
pendix |^. Below we point out some of the distinguishing 
features of the various cases. 



Case 1. With conductivities of the form in Eq. (15a), 
solving V ■ j = leads to the partial differential equation 



J ^d 2 d 2 



(l + 7f) 
£ 2 



2 



ZJ 2 



V = 0, (16) 



where 7^ is a dimensionless variable given by 7^ 



A^/S AI . Separation of variables, V(x,z) 
then yields 



dx 2 
dz 4 



k l )X = 0, 



(l + 7f) 
t 2 



d?_ 

dz 2 



Yl x k 

^J 1 



= X{x)Z{z), 

(17a) 
Z = 0. (17b) 



The solution of Eq. ( |17a| ) is X(x) = Acos(kx) + 
B sin{kx). Applying the boundary condition that no cur- 
rent enters on the sides, i.e. d x V(±L/2, z) — 0, gives 
k = mr/L with B = for even n and A = for odd n. 
When n > 0, the solution of Eq. ( |17b| ) is 

Z(z) = Ae z ^+ + Be~ z /t+ + Ce z/ «~ + Ee~ z/ t~ , (18) 

which has two length scales given by 



2 1 A? 



± 



n 

a! 



l + 7f 
i 2 

p. 



An 2 



1/2- 



(19) 



where X z = ^/E Z L 2 /7T 2 E X is a length scale occurring in 
the local limit (7* = 0). Note that if j% > (corre- 
sponding to a negative viscosity coefficient) the £'s are 



real; whereas if 7* < (corresponding to a positive vis- 
cosity coefficient) the £'s can become complex. 

Case 2. For conductivities given by Eq. ( |l5b| ), the 
corresponding differential equation is 



1 

1 2 



o 2 



0? 



V = 0. (20) 



Proceeding with separation of variables, the functions 
X n (x) are exactly the same as in Case 1, and the func- 
tions Z n {z) have the same form as in Eq. (|l8|), but the 
two length scales £± are now given by 




n 

X 2 



1 

T 2 



4» 2 7l 
£ 2 A? 



1 1/2* 



(21) 



As opposed to Case 1, this time when 7* < the £'s are 
real; and when 7! > the £'s may become complex. 

Case 3. T urni ng the example with conductivities 
given by Eq. (15c) into a partial differential equation 
yields 



2q2 



(l+7g) 
£ 2 



0? 



And separation of variables leads to 



{- 

[dx 4 

dz 2 



Yi 7 K 



f: 2 



S^ 2 

H x l 2 



0. 



7 = 0. (22) 

X = 0; (23a) 
(23b) 



We see here that this situation differs from the previous 
two in that the solution of X(x) is no longer simply sines 
and cosines, and similarly k has become nontrivial. This 
feature complicates the application of the boundary con- 
ditions and the summing over eigenfunctions necessary 
to achieve a complete solution. 

Case 4. Conductivities of the form Eq. (15d) lead to 



Q 4 + Wtftf _ J_ 5 2 



(1+7?) ,2 



P 



V = 0. (24) 



The method of solution is like that for Case 3, including 
the nontrivial values of k. 

Results. Let us now consider the results of using Pade 
approximants. We look at the behavior of V top and Vbot 
as a function of i, as in the previous section. In Cases 1 
and 2, we input a current 

hop{x) = h sin(7rx/L); hot{x) = (25) 

for the top geometry, and 

hop{x) = hot(x) = Jo [1 + sin(7ra;/L)] (26) 
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for the side geometry. These currents are chosen be- 
cause they roughly approximate the experimental inputs 
but involve the minimum number of Fourier components, 
which simplifies the calculation. If wc input the same 
currents in Cases 3 and 4, our solution will involve an 
infinite number of terms since our formula for V(x, z) is 
not in the form of a Fourier series but is a sum over more 
complicated eigenfunctions. We instead expand the cur- 
rent in terms of these eigenfunctions. In fact, we once 
again choose the input currents that minimize the num- 
ber of terms in the sum over eigenfunctions. In the top 
geometry, the current is chosen to be proportional to the 
eigenfunction corresponding to m, where K\ is the eigen- 
value which tends to ir/L in the local limit. For the side 
geometry, we choose the combination of two eigenfunc- 
tions that tends to the current of Eq. (p6f). For Case 3, 
this unfortunately means that we vary the input current 
as we vary t; at I = 0, the input current is the same as for 
Cases 1 and 2, but this smoothly evolves so that in the 
limit of t —>■ oo, the input current is It op ~ h sin(27rx / L) . 
The side geometry is affected similarly. Hence, we must 
bear in mind that the input current changes significantly 
in Case 3. In Case 4, on the contrary, the variation of 
input current with I does not appear to be as substan- 
tial. The voltages were calculated at z = 0, z = D and 
x = ±0.3L. For each case, we did two sets of calculations: 
one using the conductivities as written in Eqs. ( [l5| ) with 
A scaling as £ 2 , and the second using the scaling forms 
appro priat e f or a B ose glass; as given in Appendix ^ by 
Eqs. dC3«b — (|C3db. 



The first point to make about all of the results is that 
we never found V r i g ht/Vieft — ► 1 so long as A > 0, which 
corresponds to a negative viscosity coefficient. Hence 
none of these forms give jis the strongly nonlocal be- 
havior seen by Safajj-.pi, al. u (though perhaps they do re- 
semble other resultsaliiJ). We showed in the infinite-slab 
case (Sec. Ill), in which the conductivities corresponded 
most closely to Case 4, that we could only model Safar's 
results with positive viscosity coefficients. The results of 
this section suggest that this statement may hold for all 
of the cases. 



twinned 




3 
D/l 



FIG. 7. The results obtained for the ratio Vbot /Vtap for 
Case 2. The conductivit ies u sed are given by Eq. (15b) for 
the lower curve, and Eq. ( C3h ) for the upper curve (which has 
the appropriate scaling for a twinned sample). The parame- 



ters used are T, x — cri"' = 5, E z = a?"' — 1 



» 



I D = l, 



L = 5, Ci = 1, C2 = 5 and = 5£ 2 — note the scaling with 
£ which is present. For the twinned case, Vbot/Vtop does in 
fact tend to 1 as I — > 00. 

Secondly, we look at the difference between the Bose- 
glass scaled and unsealed results. We have already seen 
one example of this comparison in Fig. |^. Another ap- 
pears in Fig. fj], this time for Case 2. The curves are re- 
markably different considering they c ame from the same 
overall form of the conductivity, Eqs. ( |l5b| ), which shows 
how vastly different behaviors can be modeled by the 
same form, making it difficult to extract detailed infor- 
mation on the conductivity from the experimental data. 

Recall that the analysis of Huse and Majumdarlll, 
discussed in the introduction, has a symmetry under 
Vxzzx Vzxxz- If this symmetry applies in the Pade anal- 
ysis, it would correspond to a symmetry between Cases 
2 and 4. We have looked for such a symmetry. The com- 
parison is complicated by the fact that the eigenfunctions 
X n (x) are different in the two cases, as mentioned above. 
However, using conductivities given by Eq. @ led to 
qualitatively different features for Cases 2 and 4 — while 
admittedly the Case 4 input currents change in this anal- 
ysis, we do not expect the difference to affect the general 
features of the results. Thus, the Pade forms do not seem 
to share the symmetry found in hydrodynamics. Further- 
more, this lack of symmetry seems to persist even when 
we choose Pade forms with positive viscosity coefficients 
that should correspond more closely to the hydrodynamic 
case. 

Notice that in the unsealed data in Fig. [?] Vbot goes 
negative. Although the potential reverses sign, calcu- 
lations reveal that the current flow is always from left 
to right, even at the bottom of the sample. In Sec. || 
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we showed that if the conductivity is nonlocal, then ei- 
ther the real-space resistivity or conductivity is negative 
at some point. We believe that the sign reversal in the 
unsealed data is simply a consequence of this fact. It 
should be noted that voltage reversals haiia in fact been 
measured by S. Aukkaravittayapun et al.cB; however, it 
seems unlikely that nonlocal conductivity accounts for 
their results. 

Before concluding this section, we note that it is possi- 
ble to extend the Pade forms considered above to higher 
levels of approximation. Suppose one simply adds a sec- 
ond Pade term to the first one, for instance, 



1 + klt\ 1 



1,2/2 1 



(27) 



There are now two length scales, and one can choose the 
parameters so that the conductivity has a positive viscos- 
ity and maintain the property er^^O) > o xx {k x — ► oo). It 
is still possible to solve the equation V • j = in a man- 
ner similar to the examples discussed above. However, 
one would need to differentiate the integro-differential 
equation four times instead of twice in order to eliminate 
the two integral terms. As a result, one ends up with 
a sixth-order equation and a lot more algebra. We have 
not pursued this avenue. 



V. SURFACE CONSIDERATIONS 

We begin this section on surface considerations with a 
comparison of the voltage and current distributions found 
using the hydrodynamic approach with those found us- 
ing the Pade approach. Let us consider Case 2 from the 
previous section making the following parameter choices 



Ex = °l - 
A z = 



- v r 2 7 



(28a) 
(28b) 



so that the small-k expansion of the conductivity is 

a xx (k) = ol + V kl - r]t 2 ki + 0(£ A k 6 z ). (29) 

Note that the £ — > limit coincides with the truncation 
used in the hydrodynamic approach. The fourth-order 
equation derived for Case 2, Eq. (EG), becomes 



d 2 z dl-a° x d 2 x -X z d*]V = 0. (30) 



If the terms proportional to £ 2 are dropped, the equa- 
tion is identical to the one studied by Huse and Ma- 
jumdarli3c3. However, notice the small parameter, £ 2 , 
multiplies the highest derivative, d z V. This is the clas- 
sic scenario for the development of a boundary layer, a 
small region in which V varies-quite rapidly and in which 
£ 2 d z V is not negligibly small O Recall that our analysis 
of Case 2 produced two length scales given by Eq. (21). 
In the small-£ limit, £ + is a small length (proportional 



to £); while £_ is the length scale arising in the hydro- 
dynamic analysis. It is tempting to conclude that one 
has discovered the length scale associated with Huse and 
Majumdar's surface currents, but a more careful analysis 
is in order. 

Following the steps outlined in the Appendix |b] we 
can calculate V(x, z) for the top geometry with an in- 
put current 1\ smijtx j U). ^From it we calculate j x (x, z), 
the current density in the x-direction. Figure |^ shows 
j x (0,z) as a function of z for several values of £. As £ 
decreases, the current becomes increasingly confined to 
the surface z = D, in other words, we obtain surface cur- 
rents. However, they differ from those found by Huse and 
MajumdarcJ. For the same input current the hydrody- 
namic analysis predicts that the ratio of current carried 
in the surface to that in the bulk is 



surface current 
bulk current 



(31) 



whereas Fig. || seems to show that all of the current in 
our solution is carried in the surface in the limit £ — > 0. 
We confirm this result by noting that in the bulk j x ~ £ 
in the small- £ limit, implying that we have no bulk and 
all surface current as £ — > 0. 




z/D 



FIG. 8. The current density jx(0,z) for different val ues o f 
£ corresponding to the Case 2 Pade conductivities, Eqs. (15t) 
and (Eq), with parameters: <r° = 5, T] = 5, = 1, D = 1, 
L = 5 and input current I(x) = Iisin(nx/L). The dashed 
curve is the local curve (£ — > oo), and the other curves are for 
t — 0.4, and I — 0.067. The areas under the curves are equal, 
but as £ becomes smaller there is more current density at the 
surface. 
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We have examined the I — > limit not because we 
believe it to model the real conductivity but in order to 
make contact with and perhaps better understand the 
hydrodynamic theory. Nevertheless, the outcome — that 
much of the current is confined to the surface — has been 
suggested in other contexts. If a significant fraction of 
the current flows in the surface, one might expect to find 
nonlinear behavior down to very low currents, since the 
current density near the surface would vary quite rapidly, 
leading to enormous tearing forces on the vortices. Non- 
linear IV characteristics, have indeed been observed in 
both the top geometryc3 and in c-axis resistivity mea- 
surementsE2l. However, our work is concerned only with 
the linear regime, and so we do not discuss this further. 

The fact that the small-^ limit of Eq. ( ^o|) is identi- 
cal to the corresponding hydrodynamic equation and yet 
the two approaches predict differing amounts of surface 
current suggests that it is in the treatment at the bound- 
ary (or in the effect of the boundary layer) that the two 
approaches differ. In the Pade calculation we considered 
the effect of the surface only through the boundary con- 
ditions on the current; we neglected any effect the sur- 
face might have on the nonlocal conductivity itself. One 
expects some surface dependence since the conductivity 
is determined by the superconducting order parameter, 
which in turn depends on the boundaries. Within the 
limits of their calculation, Blum and MooreEJ gave an 
explicit expression for the conductivity in the presence 
of a surface. In addition to the usual bulk conductivity 
their analysis yielded a term corresponding to the im- 
age of the bulk conductivity as well as "cross terms." In 
order to proceed with their voltage-distribution calcula- 
tion analytically, the cross terms were dropped with an 
argument suggesting their effect was small. Subsequent 
workE!J has shown that the effect while small propagates 
farther into the bulk than was suggested in that work, 
indicating again the importance of treating the surface 
effects properly. 

It may turn out that the hydrodynamic approach actu- 
ally incorporates some of these surface effects. We have 
found that the small-£ limit of the Pade approach with 
some surface effects duplicates the hydrodynamic results. 
In the previous analysis, we used the conductivity in 
Eq. (|l5bp which corresponds to 



a xx (r,r') = 



yG(z,z') 



S(x 



with 



G(z,z') 



J_ -\z-z'\/t 

21 



Now, we modify this conductivity so that 



G(z,z') = 



cosh 



D-\z-z'\ 



cosh 



D-z~z' 



21 sinh[L>/^ 



x') 
(32) 

(33) 
(34) 



which corresponds to including a series of image terms 
such that the nonlocal conductivity satisfies the bound- 
ary condition of its derivative vanishing on the two sur- 
faces z = and z — D. Some motivation for using this 
choice might come from the—boundary condition the or- 
der parameter itself satisfiestll; however, the justification 
here is that the results match those of the hydrodynamic 
approach. 

Using this form, one can still solve the V • j = 
equation in the same manner as for the conductivity of 
Eqs. (|l5|). In fact, the partial differential equation one 
obtains as an intermediate step is the same as that ob- 
tained in Case 2 with a translationally invariant conduc- 
tivity, i.e. Eq. (20). But this is to be expected since the 
difference lies in the boundary and not in the bulk. Ap- 
plying the boundary conditions results in the following 
voltage distribution 



V(x,z) = 



tit- (e-^)cosh(^ 



^zp {e+ 



£M sinh 



(35) 



which agrees with the hydrodynamic result in the I — > 
limit. The corresponding current distribution j x (x,z) 



jx(x,z) 



hi cob(^)££_ (V - <J Q x e + <e) cosh(^) 



+{£- 



£2) silm (_D j 



(36) 



is plotted in Fig. Q. It can be seen that once again we 
have surface currents in the limit I — > 0; however, this 
time only some of the current flows in the surface, with 
the rest still flowing in the bulk of the sample. In fact, 
the ratio of surface current to bulk current is identical to 
that from hydrodynamics. 
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FIG. 9. The current density j x (0, z) for different values of 
I corresponding to nontranslationally invariant conductivities 
of the form given by Eq. ( p^ ) with the same parameters and 
input current as in Fig. ^. The upper dashed curve is the local 
curve (I — * oo), and the curves with »'s are for I — 0.4 and 
I — 0.067. The lower solid curve is the Huse and Majumdar 
result (jx(0,z) ~ 0.48 cosh(0. 8152:) for the parameters used 
here) , which corresponds to the £ — » limit of our results. 



the importance of the surface in the problem of nonlocal 
conductivity and suggests that even getting the correct 
large-fc behavior of the bulk conductivity may be insuf- 
ficient if one has neglected surface effects. Fortunately, 
we have at least one example in which including the sur- 
face effects did not destroy the solubility of the Pade 
approach. 



This result is interesting in that it suggests that Huse 
and Majumdar's boundary condition may be related to 
surface effects in the conductivity. It is surprising that 
these two approaches: hydrodynamics with its transla- 
tionally invariant conductivity and its slightly unusual 
boundary condition (associated with the discontinuity in 
the derivative of the electric field at the surface) and the 
Pade approach with a non-translationally invariant con- 
ductivity and ordinary boundary condition, can produce 
the same outcome. 
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VI. CONCLUSIONS 



APPENDIX A: ANALYSIS OF THE 
INFINITE-SLAB GEOMETRY 



It is clear from the above work that the problem of 
nonlocal conductivity is not completely understood; how- 
ever, we have made progress in developing a calculational 
scheme which appears to be more robust than what was 
available previously. Firstly, with Pade approximations, 
we can solve for voltage distributions having either pos- 
itive or negative viscosities; whereas, hydrodynamics is 
limited to positive viscosities. Moreover, these approxi- 
mations give the correct behavior at both high- and low- 
k and so may model the surface effects better than the 
hydrodynamic forms. Because measurements in the stan- 
dard flux-transformer geometry are taken at the surface, 
an approach that deals inadequately with the surface has 
got to be considered suspect. 

Since we now have the ability to investigate both pos- 
itive and negative viscosities, we can make deductions 
about which sign Yields the strongly nonlocal behavior 
seen by Safar et al. B We have found that modeling these 
effects seems to require positive viscosities, particularly 
so when the nonlocal behavior in the x direction in a zz 
is considered. As discussed in Sec. |l[ there is mounting 
evidence that cr zz {k z ) may have a negative viscosity coef- 
ficient even at low temperatures, implying that the non- 
local behavior seen is not due to the dependence of a zz 
on k z . (There are, however, other current configurations 
in whictuthis is the only dependence of the conductivity 
probed £3) 

Finally we have discussed how the conductivity itself 
may be affected by surfaces and have shown that remov- 
ing the translational invariance of the conductivity does 
not necessarily complicate the voltage-distribution calcu- 
lation. In particular, we have shown that the Huse and 
Majumdar solution corresponds to the limit of a Pade- 
type solution in which the conductivity lacks translation 
invariance. This surprising outcome emphasizes again 



In this appendix we consider the length scales and 
surface effects that characterize Vr(x,z) and Vs(x, z), 
the voltage distributions in the infinite-slab geometry, by 
considering the pole structure of the integrals given in 
Eqs. ( p"2"| ) and (O) fo r a number of choices of a 7 , z (k x ) . The 
Vs integral, Eq7(|l3[) , is done by summing the residues as- 
sociated with the zeros of cosh(KD/2), given by n(k x ) = 
±z7t(2to + 1)/D with m = 0, 1, . . .; the Vr integral is 
rather similar. For the local problem (jj zz (k x ) = T, z ) 
these poles are evenly spaced along the imaginary axis. 
In fact, with <5-function-distributed input currents 

J T = J [ 5(x - 2L C ) - 5{x + 2L C )] (Ala) 
Js = -J d(x + 2L C ), (Alb) 

the resulting expressions can be resummed to give 



V T , s {x,z)=V In 



cosh 



2L C + x 



cos 



D 



cosh ± cos — 



(A2) 



where V = — Jo/2n\/'S x E 2 and \ x — ^Y, X D 2 /tt 2 T, z and 
where the upper signs in the denominator correspond to 
Vt{x, z) which is odd about x — and the lower signs 
correspond to Vs(x, z) which is odd about z — D/2EJ. 

If we consider a nonlocal part to have the hydrody- 
namic form a zz (k x ) = cr z + rjk 2 ,, the poles are located 
at 



k x (m) 



±i(2m+ 1) 



A2+77(2m+l) 2 /<7° 



1/2 ' 



(A3) 
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where \ x = ^JYj x D 2 / '7r 2 er°. Note that the smallest pole 
(to = 0) is shifted to smaller k x compared to the local sit- 
uation (j] = 0), implying a longer length scale. There is 
only a significant shift if the viscous length scale \/ri/a z 
becomes comparable to the "local" length scale A^, which 
depends on the sample thickness D. As 77 increases fur- 
ther the viscous length dominates; we were able to resum 
the series in this limit, finding 



V s (x,z) 



- J (2z - D) 



■ exp 




(A4) 



It would be interesting to probe the spatial dependence 
of Vs experimentally and extract its length scale. How- 
ever, such a measurement would be difficult as it would 
require a sample long enough to accommodate several 
leads, and the voltages far from the current may become 
too small to be meaningful. Returning to the pole struc- 
ture in the hydrodynamic case, another point to notice 
is that they accumulate at a finite value ±iy/a^Jrj. As a 
result Vs(x,z) no longer diverges logarithmically at the 
contact points (—2L C ,D) and (— 2L C ,0). Given the 6- 
function input currents, this divergence is physical, and 
the failure of the hydrodynamic form to reproduce it is 
an example of how the incorrect large-k can affect the 
potential especially at the surface. 

Next, let us consider a zz (k x ) to have a Pade form 



s (k x ) 



1 + k 2 £ 2 ' 



(A5) 



The small-fc^ behavior of this expression is similar to the 
hydrodynamic example if < 0. However, it has twice 
as many poles since n{k x ) = ±i7t(2to + l)/D has twice 
as many solutions as in the hydrodynamic case. These 
poles break into two sets. For large to one set is evenly 
spaced and mimics the behavior in the local problem in- 
cluding the logarithmic divergence at the leads, while 
the other set accumulates at the value ±i-v/(T + 7* )/£ 2 
(where 7J = Af/£ z ) and mimics the hydrodynamic be- 
havior in the bulk. The Pade form has clear advantages 
over the hydrodynamic form, but even it is not quite 
right since the large- k x limit of a zz (k x ) should be less 
than a zz (0). To achieve that one needs something like 



1 



71 



72 



1 + k 2 J 2 (1 + k 2 J 2 ) 2 



(A6) 



with 71/2 < 72 < 71, which leads to three sets of poles — 
one like the local case and two lik e the hydrodynamic 
case. With the Pade form, Eq. (A5), one can also inves- 
tigate the consequences of "negative" viscosities when 
Af > 0. As Af increases, the i-axis length scale de- 
creases and eventually a point is reached at which the 
poles move off the purely imaginary axis and some oscil- 
latory behavior is superimposed on the exponential de- 
cay. 



APPENDIX B: SOLUTION OF THE 
INTEGRO-DIFFERENTIAL EQUATION 
ASSOCIATED WITH THE PADE FORM 

In this appendix we outline the solution of Case 2. For 



conductivities of the form given by Eq. (15b), the steady- 



state continuity equation V • j = takes the form 



V x d z x V{x,z) + Y, z d z z V(x,z) 



A z 
' 21 



-I* 



z '\ /t dlV{x,z')dz' 



0. 



(Bl) 



z'\/£}/2£ is the Fourier transform of 



where exp{— \z 
(l + k 2 J 2 )-\ 

Le t u s now apply the following trick. Differentiate 
Eq. (Bl) twice with respect to z, which leads to 



T, x d 2 z d x V{x, z) + Y, z d%V(x, z) 



A 

T 2 



x %V(x,z) 




e- lz - z ' l/e d 2 x V(x, z')dz' \ = 0, (B2) 



where we have exploited the relation 
d 2 



dz 2 



\z-z'\/t\ _ 



-\z-z'\/e 



(B3) 



familiar from the solution of Schrodinger's equa tion with 
a (^-function potential. We can combine Eqs. (Bl) and 
(B2) to eliminate the integral term, giving 



v z dt + z x d 2 z d 2 x 



(Sx + A|) 2 _ E, 2 
I 2 £ 2 



V = 0, (B4) 



which was given in the main body of the paper as 
Eq. ©. 

Separating variables, V(x,z) — X[x) Z(z), and ap- 
plying the boundary condition j x (±L/2, z) = yields 
X(x) = Acos(2mrx/L) or X(x) = Bsm[(2n + l)irx/L] 
for n = 0, 1, . . .. When n > 0, the corresponding Z(z) is 
given by 




(B5) 



where £± = £± (n) is given by Eq. (|T|) . For some of the 
algebra that follows it is convenient to choose modes that 
are symmetric about z = D/2. 

It might appear that we need to apply four boundary 
conditions to determine the constants in Z n (z), but ac- 
tually two of the constants are found by inserting the 
solut ion into the original integro-differential equation, 
Eq. (Bl). This step yields the following two conditions 
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Pn 
Rn 



d-(£ 2 ~e + ) ^cosh^+^sinh^) 



£+cosh(^)+£sinh(^; 



fcosh^ 



■ e_ S inh(^; 



£ + sinh(^; 



, (B6a) 



(B6b) 



The two remaining constants are fixed by the boundary 
conditions on j z , namely 



X z d z V(x,D)=I top (x), 
-Z z d z V(x,0)=I bot (x). 



(B7) 



They turn out to be fairly complicated functions of the 
parameters £, E^, E z , 7^, n and L, in addition to the 
Fourier components of Itop{x) and I bo t{x). Note that 
each Z n (z) is only a function of the corresponding Fourier 
component of the current. The n — part of the solution 
requires separate consideration; however, it is straightfor- 
ward and no details are provided here. 



APPENDIX C: BOSE GLASS SCALING FORMS 
FOR CONDUCTIVITIES 



In this section we consider scaling forms for the con- 
ductivities which might be expected to hold in the pres- 
ence of correlated disorder {e.g. columnar defects or twin 
boundaries). With such disorder, there is thought to 
be a second-order phase transition between the low- 
temperature Bose-glass phase and the high-temperature 
phase .consisting of an entangled liquid of delocalizcd flux 
linesafl. Near the transition the characteristic length 
scales £± (within the ab planes) and £u (along the c-axis) 

and the characteristic time scale r ~ £\_ diverge. Nelson 
and RadzihovskyEa used the scaling of the free-energy 
density, / ~ l/£n£j_, and the vector potential, Am ~ 
and A± ~ l/£±_ (from gauge invariance), and the rela- 
tions J = df /dA and E = —dA/dt to suggest that the 
conductivities scale as follow 
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f 2 



4 



(C2) 



where the superscript (s) refers to the superconducting 
contribution, we also include in T, x and E z normal contri- 
butions Ox and o z n ^ that are not affected by the scaling. 
We thus obtain the following Bose-glass scaling forms: 
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It should be noted that although the constants C\ and C 2 
have no explicit dependence upon the length scales, they 
will be temperature-dependent. However, compared to 
the temperature dependence of the length scales £± and 
£n near the transition, which go as powers of \T — Tbg\ 
(where Tbg is the transition temperature), it is a weak 
dependence. In the forms above the coefficient of the 
Pade term is assumed to be positive, and therefore the 
viscosity coefficient is negative. We can also write forms 
that have positive viscosity coefficients; for instance, Case 
4 would be 



5 (k) 
: (k) 



= <>) 



» 



C 2 £i, 
2C X £\ 



+ k^£\ 



(C4) 



which is the form used to generate Fig. |^. 



Studies of this transitionQ'L 2 ] have suggested that £11 ~ £\ 
and z' — 6.0 ± 0.5; we are going to use z' = 6. 



In the Pade-form conductivities (Eqs. 15), we have a 
length scale £, which is a c-axis length scale £h in Cases 1 
and 2 and is an a6-plane length scale £±_ in Cases 3 and 
4. Recall the conductivities in Case 1 are o xx {k) = E T 
and o zz (k) = E 2 + A*/(l + k 2 £ 2 ). What we want to do 
here is determine the dependence of the constants E and 
A upon £. Since in this case the length scale is £\\ , we use 
£11 ~ £j_ to eliminate the dependence upon £± and arrive 
at 
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